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Using the Hartle-Hawking no-boundary proposal for the wave function of the universe, we can 
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wave function and probability of a universe with a black hole are calculated at the WKB level. The 
probability of a black hole creation is the exponential of one quarter of the sum of areas of the black 
hole and cosmological horizons. One quarter of this sum is the total entropy of the universe. We 
show that these arguments apply to all kinds of black holes in the de Sitter space background. 
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I. Introduction 



Hawking's theory of No-Boundary Universe has for the first time led to a self-contained cosmology. 
Now in principle, one can predict everything in the universe solely from physical laws and the long- 
standing first cause problem has been dispelled. In quantum cosmology one of the most challenging 
problems is the existence of primordial black holes. 

It is well known that a black hole can be formed in two ways, the first being the gravitational 
collapse of a massive star. If the mass of a star exceeds about twice that of the Sun, a black hole 
will be its ultimate corpse. The second way originates from the fluctuation of matter distribution 
in the early universe. In the big bang model, the matter content can be classically described [1][2], 
while in the inflationary universe the matter content is attributed to the quantum fluctuation of the 
Higgs scalar [3]. 

The discovery of Hawking radiation of black holes has been the most important event of grav- 
itational physics for several decades. However, the life of the first kind of black hole is very much 
longer than the age of the universe. The only hope of confirming Hawking radiation by observation 
is through primordial black hole hunting. 

Strictly speaking, black holes formed through either way mentioned above can hardly be regarded 
as primordial. A true primordial black hole should be created at the moment of the birth of the 
universe. Over the last decade there have been several attempts to deal with this problem, however 
their results are not conclusive [4] [5] . 

It is believed that the very early universe is approximately described by a de Sitter metric. In 
quantum cosmology, at the Planckian era, the universe was created from a S 4 space through a 
quantum transition. Therefore, to study the problem of primordial black hole creation in the de 
Sitter spacetime background is of twofold interest, for cosmology and for black hole physics. 

There have been many studies recently on quantum creation of charged or neutral black hole pairs 
in the de Sitter spacetime background [6] [7] [8] [9] [10] [II]. The case of a single primordial black hole 
is the topic of this paper. Sect. II is devoted to the theory of constrained gravitational instanton. 
Sect. Ill considers both the neutral and nonrotating black hole, i.e. the Schwarzschild black hole 
and the charged but nonrotating black hole, i.e. the Reissner-Nordstrom black hole. Sect. IV is 
devoted to the rotating but neutral black hole case, i.e. the Kerr black hole. Sect. V investigates 
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the rotating and charged case, i.e. the Newman black hole. By the no-hair theorem, all these kinds 
of black holes have exhausted the stationary vacuum or electrovac cases. Therefore, the problem of 
quantum creation of a single black hole in quantum cosmology is completely resolved. Sect. VI is a 
discussion. 



II. The constrained gravitational instantons 

In the No-Boundary Universe the wave function of the universe is given by [12] 

V{h ij ,<f>)= d[ff M „]d[0]exp(-J([ff M ,,,$), (1) 
Jc 

where the path integral is over class C of compact Euclidean 4-metrics and matter field configura- 
tions, which agree with the given 3-metrics hij of the only boundary and matter configuration <p on 
it. Here I means the Euclidean action. 

The Euclidean action for the gravitational part for a smooth spacetime manifold M with bound- 
ary dM is 

/=-— / d i xq l ' 2 {R-2k)-— ( d 3 xh}l 2 K, (2) 

where A is the cosmological constant, R is the scalar curvature, K is the trace of the second funda- 
mental form of the boundary, g and h are the determinants of g^ and hij respectively. 

The dominant contribution to the path integral comes from some classical solutions of the field 
equations, which are the saddle points of the path integral. 

The probability of the Lorentzian trajectory emanating from the 3-surface £ with the matter 
field 6 on it can be written as 



P = = / d[g^]d[4>] exp(-/([ ff ^,0]), 
Jc 



(3) 



where class C is all no-boundary compact Euclidean 4-metrics and matter field configurations which 
agree with the given 3- metric h^ and matter field (f> on S. 

Here, we do not restrict class C to contain regular metrics only, since the derivation from Eq. 
(1) to Eq. (3) has already led to some jump discontinuities in the extrinsic curvature at S. 

The main contribution to the path integral in Eq. (3) is due to the stationary action 4-metric, 
which meets all requirements on the 3-surface X and other restrictions. At the WKB level, the 
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exponential of the negative of the stationary action is the probability of the corresponding Lorcntzian 
trajectory. 

From the above viewpoint, an extension of class C to include metrics with some mild singularities 
is essential. Indeed, in some sense, the set of all regular metrics is not complete, since for many 
cases, under the usual regularity conditions and the requirements at the equator E, there may not 
exist any stationary action metric, i.e. a gravitational instanton. It is not clear, how large class C 
should be. A necessary condition for a metric to be a member it that its scalar curvature should be 
well-defined mathematically. It is reasonable to include jump discontinuities of extrinsic curvature 
and their degenerate cases, that is the conical or pancake singularities. For this kind of singularity, 
the quantity g x ^ 2 R can be interpreted as a distribution- valued density [13]. 

Although the regularity conditions on the 4-metrics and the requirements from the equator 
E sometimes are so strong that no gravitational instanton exists, one can still hopefully find a 
stationary action nonregular solution with some mild singularities within class C, which can be 
called the constrained gravitational instanton. Here the manifold of the instanton is constrained by 
the equator E. 

It has been proven [13] a stationary action regular solution keeps its status under the extension 
of class C . However, if a stationary action regular solution cannot be found, then it can probably 
be expected with some singularities in class C . For a model with S 1 x S 2 topology under the 
minisuperspace ansatz 

ds 2 = a 2 (r)dr 2 + b 2 (r)dr 2 + c 2 (r)rffi 2 , (4) 

where z is periodic with period 2ir, and dfl?, represents the metric of a unit 2-sphcre, the solution 
satisfies the usual Einstein field equation except for the singularities at the final r — rj and initial 
r = Ti surfaces. One can rephrase this by saying that the solution obeys the generalized Einstein 
equation in the whole manifold. Since this result is derived from first principles, one should not feel 
upset about this situation. 

Except for the interpretation of probability the above arguments can also be applied to the 
Lorentzian regime with a purely imaginary phase. The dominating contribution is again due to the 
stationary action trajectories. However, in most cases, the restrictions are not too strong, and one 
can find a regular metric satisfying the usual Einstein field equation. 
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At the transition surface S, it is assumed that along neither of the sides of S does a singular 
matter distribution exist. It follows from the Einstein equation that the fundamental form Kij at 
£ should vanish, 

KH = 0. (5) 

This condition cannot apply to the mild singularities at £ if there is any, since the usual Einstein 
equation does not hold there. 

The singularity problem associated with a gravitational instanton is not always disturbing; in 
fact it can be beneficial. If the restrictions are weak enough to allow a regular instanton, then the 
Lorentzian evolution originating from it must be most probable one. Therefore, in order to find the 
most probable Lorentzian evolution, one needs only find a regular instanton, which then identifies 
the 3-metric and matter field on S. 

In general, the wave packet of a wave function of the universe represents an ensemble of classical 
trajectories. Under our scheme, the most probable trajectory associated with an instanton can be 
singled out [14]. Thus, quantum cosmology obtains its complete power of prediction. It means there 
is no more degree of freedom left as long as the model is well-defined. 

On the other hand, the more severe the restrictions are, the larger the stationary action is, and 
therefore, the less probable its corresponding Lorentzian evolution. This is the situation with a 
constrained instanton. Wc shall see this in the case of a primordial black hole. 

If there is no black hole in the universe, then one can get a regular instanton S 4 . If there is, then 
the restrictions are strong enough to forbid regular solution. Therefore, the probability of a universe 
without a black hole is always greater than one with a black hole. Our calculation will support this. 

There has been some progress in this direction. However, nearly all scenarios studied are asso- 
ciated with pair creation of black holes [6] [7] [8] [9] [10] [11]. The main reason for this is that, people 
consider our universe to have been created by a quantum transition from a gravitational instanton. 
There does not exist any gravitational instanton which provides the seed for the creation of a single 
black hole in the de Sitter background. 

As we mentioned above, in quantum cosmology one uses a Lorentzian metric to join a Euclidean 
metric, both being sectors of a complex manifold. However, there exist very few complex manifolds 
satisfying the Einstein equation with both a Euclidean and a Lorentzian sectors [15]. One may 
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appeal to some approximately Euclidean or Lorcntzian sectors, but only at the price of losing some 
of the beauty of the theory. In the extended framework the requirement becomes quite loose. The 
situation of black hole creation we are going to investigate is the best illustration. 



III. The spherically symmetric black hole 

Let us begin with a quantum spherically symmetric vacuum or electrovac model with a positive 
cosmological constant A. The cosmological 

constant may be effective due to the Planckian inflation in the Hawking massive scalar model 
[16]. At the semiclassical level the evolution of the universe is described by its classical solutions. 
The Schwarzschild-de Sitter spacetime with mass parameter m and zero charge Q is the unique 
spherically symmetric vacuum solution to the Einstein equation with a cosmological constant A. 
The Reissner-Nordstrom-de Sitter spacetime, with mass parameter m, nonzero charge Q and a 
cosmological constant A, is the only spherically symmetric electrovac solution to the Einstein and 
Maxwell equations. Its Euclidean metric can be written as [17] 

2m Q 2 Ar 2 \ _,_ 2 ( 2m Q 2 Ar 2 \ _1 2 2 2 . 2 2 



&=[!- — + ^-—\df+ll- — + ^- — ) ^+rW+sinW). (6) 



We can set 

2m Q 2 Ar 
r r 2 3 
For convenience one can make a factorization 



V. = l-™ + *t-^. (7) 



V s = -7^( r - r o)(r - ri)(r - r 2 )(r - r 3 ), (8) 

where ^Oj r ij r 2, ^3 are in ascending order. r 2 and r% are the black hole and cosmological horizons, 
where conical singularities may occur, tq is negative. If the black hole is neutral, then r\ can be set 
to zero, and there are essentially three roots left. 
The gauge field is 

F = - l -QdrAdr (9) 

for an electrically charged solution, and 

F = Qsin9de Adct) (10) 
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for a magnetically charged solution. We shall not consider dyonic solutions. 
The roots satisfy the following relations 



£^ = -f' (12) 



i>j 



A 



E6m , „. 

nr r k = -— (13) 

i>j>k 

and 

n^=-¥- (14) 

The black hole and cosmological surface gravities K2 and ^3 are [13] 

«2 = ^'(ra)! = A (r 2 - r )(r 2 - n)(r 3 - r 2 ), (15) 



«3 = ^'faOl = ^2 (r 3 - r„)(r 3 - n)(r 3 - r 2 ). (16) 



1, A_ 

The requirement of vanishing second fundamental form at £ minus the two conical singularities 
at the two horizons implies that the transition can only occur at two sections of constant values 
of imaginary time r glued at the two horizons. The 3-surface S has topology S 2 x S 1 . To form 
a constrained gravitational instanton, one can have two cuts at r = const s. between r = r 2 and 
r = r 3 . Then the /2-fold cover turns the (r — r) plane into a cone with a deficit angle 2n(l — / 2 ) at 
the black hole horizon. In a similar way one can have an /3-fold cover at the cosmological horizon. 
Both / 2 and / 3 can take any pair of real numbers with the relation 

/ 2 /?2 = / 3 /?3, (17) 

where /3 2 = 27tk 2 ~ 1 and /?3 = 27TKg 1 . If / 2 or / 3 is different from 1, then the cone at the black hole 
or cosmological horizon will have an extra contribution to the action of the instanton. After the 
transition to Lorentzian spacetime, the conical singularities will only affect the real part of the phase 
of the wave function, i.e. the probability of the creation of the black hole. 
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Since the integral of K with respect to the 3-area in the boundary term of the action (2) is the 
area increase rate along its normal, then the extra contribution due to the conical singularities can 
be considered as the degenerate form shown below 

^deficit = ~ ■ 47rr| • 2tt(1 - / 2 ), (18) 

h,deficit = ~ ■ 4tit§ • 2tt(1 - / 3 ). (19) 
The action due to the volume is 

h = -^(rl - 4) ± Mf-tf - (20) 

where + is for the magnetic case and — is for the electric case. This term disappears for the neutral 
case. 

In the neutral case, the boundary date on the 3-surface X will be hij. In the magnetic case, the 
boundary date is hij and Ai. The vector potential in turn determines the magnetic charge, since it 
can be obtained by the magnetic flux, or the integral of the gauge field F over the S 2 space sector. 
It is more convenient to choose a gauge potential 

A = Q(l-cos6)d<j> (21) 

to evaluate the flux. 

In the electric case, the boundary date is hij and the momentum u [11], which is canonically 
conjugate to the electric charge and defined by 



(22) 



where the integral is around the S 1 direction. The most convenient choice of the gauge potential for 
the calculation is 

A=- l -9- T dr. (23) 

The wave function for the equator is the exponential of half the negative of the action. For the 
neutral and magnetic cases, one obtains the wave function ^(/ly) and For the electric 
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case, what one obtains this way is ^(w, h%j) instead of ^(Q,hij). One can get the wave function 
*(Q, hij) for a given electric charge through the Fourier transformation [10] [11] 

1 f°° 

*(Q.M = ^-/ dwe*"«*(o;,/iy). (24) 
This Fourier transformation is equivalent to a multiplication of an extra factor 

ex P (-f^Wf-'i 1 )) ( 25 ) 
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to the wave function. This makes the probabilities for magnetic and electric cases equal, and thus 
recovers the duality between the magnetic and electric black holes [11]. 

Finally, using the relations (17) and (11)-(14), one obtains the probability for a spherically 
symmetric black hole creation 

P s ^exp(Tr(r 2 2 +rD). (26) 

This is the exponential of one quarter of the sum of the black hole and cosmological horizon areas, 
or the total entropy of the universe. 

The most remarkable fact is that the result is independent of our choice of ji or f^. It means the 
manifold has a stationary action, therefore it can be qualified as a constrained gravitational instanton, 
and it can be used for the WKB approximation to the wave function. The same phenomenon will 
occur to the Kerr-Newman case as one will see later. 

For the cases of the nonsingular, charged or neutral, spherically symmetric instantons and the 
associated black hole creations [6] [7] [8] [9] [10] [11], all these instantons lead to the creation of pairs of 
black holes. For these cases one can avoid the conical singularities by choosing fi = = 1, since 
the two surface gravities are identical. However, their results are the special cases of our general 
formula (26), recalling that the degenerate horizon should be counted twice. 

The wave function for the spherically symmetric black hole can also be found [5]. 

When m = and Q = 0, it is reduced to the de Sitter case 

P QKiexp (^.) (27) 



^ A , 

and when Q = and r 2 = r 3 , it is reduced to the Nariai case 
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The formula (26) interposes the above values for the two extreme cases of neutral black holes. 

The probability is a decreasing function with respect to parameter m and \Q\. So the de Sitter 
universe is the most probable one for the Planckian era in quantum cosmology, as is expected. 



IV. The Kerr-de Sitter black hole 



Now let us discuss the creation of a rotating black hole in the de Sitter space background. The 
Lorentzian metric of the black hole spacetime is [17] 

ds 2 = p 2 {A- 1 dr 2 + A~ g 1 d9 2 ) + p- 2 ~r 2 A B sin 2 9(adt - (r 2 + a 2 )dcp) 2 - p- 2 ~- 2 A r {dt -asm 2 9d<f)) 2 , 

(29) 

where 

p 2 = r 2 + a 2 cos 2 9, (30) 

A r = (r 2 + a 2 )(l - Ar^ 1 ) - 2mr + Q 2 + P 2 , (31) 

A e = 1 + Aa 2 3 _1 cos 2 6», (32) 

S^l + Aa^ 1 (33) 

and m, a, Q and P are constants, m and ma representing mass and angular momentum. Q and P 
are electric and magnetic charges. 
One can factorize A r as follows 

A r = ~(r - r )(r - n)(r - r 2 )(r - r 3 ), (34) 

where the roots r , r 1; r 2 and r 3 are in ascending order, r 2 and r 3 are the black hole and cosmological 
horizons. The roots satisfy the following relations: 



5> = ' ( 35 ) 



^nrj = -j+a 2 , (36) 

i>j 
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E6m 
nrj-rfc = - — , (37) 

n-.=- 3( ° 2+ f +p2) . OS, 

In this section we shall concentrate on the neutral case with Q = P = 0. The Newman case with 
nonzero electric or magnetic charge will be differed to the next section. 

The probability of the Kerr black hole creation, at the WKB level, is the exponential of the 
negative half of its corresponding constrained gravitational instanton. The only instanton which can 
be used to join the Lorentzian sector at the quantum transition is the complex spacetime obtained 
from the Lorentzian metric by a substitution t — ► — it only. However, for convenience of calculation, 
we can let a to be imaginary, and then the complex metric becomes Euclidean. After we get the 
probability for the imaginary a value, then we can analytically continue back to real a to obtain the 
required probability. 

In order to form a constrained gravitational instanton, one can do the similar cutting, folding 
and covering at both the black hole and cosmological horizons with / 2 and f 3 satisfying relation (17) 
as in the nonrotating case. We shall freely switch back and forth between the real and imaginary 
values of a in the following calculation to facilitate our interpretation. 

For the Kerr case, the topology of 3-surface £ is S 2 x S 1 . Their horizon areas are 

A 2 = Mr 2 2 + a 2 )-- 1 , (39) 

A 3 = 4Tr(r 2 3 + a 2 )-- 1 . (40) 

The black hole and cosmological surface gravities are 

A(r 2 -r )(r 2 - n)(r 3 - r 2 ) , . 

K2 = c-l 2 , 2\ ' ( 41 ) 

6^(r| + a 2 ) 



A(r 3 - r )(r 3 -n)( r 3 - r 2 ) 
6E(r 2 + a 2 ) 

The actions due to the conical singularities are 



<OJV3 '1A'3 '27 

= ^2T^ ' ( 42 ) 



- ^r 2 2 +« 2 )(l-/ 2 ) 

J-2,deficit — = , (40 ) 

h.deficit = — 2- 1 — ■ (44) 
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The action due to the volume is 



where fa is defined as before. 

If one naively takes the exponential of the negative of half the total action (after the analytic 
continuation by the replacement of b by a), then the wave function for the creation moment of a 
black hole with parameter m and a will not be obtained. The physical reason is that what one 
can observe is only the angular differentiation, or the relative rotation of the two horizons. This 
situation is similar to the case of a Kerr black hole in the asymptotically flat background. There one 
can only measure the rotation of the black hole horizon from the spatial infinity. To find the wave 
function for the given mass and angular momentum one has to make the Fourier transformation 

1 f°° -2 

^(m,a,hij) = — / dSe lSJ " V(m,6,hij), (46) 

where 5 is the relative rotation angle for the time period / 2 /?2 , which is canonically conjugate to the 
angular momentum J = ma; and the factor S -2 is due to the time rescaling. The angle difference 
5 can be evaluated 

5= / dT(n 2 -n 3 ), (47) 

Jo 

where the angular velocities at the two horizons are 

* = W^- < 48 > 

and 

fi3 = ^) 
The Fourier transformation is equivalent to adding an extra term into the action for the con- 
strained instanton, and then the total action becomes 

7= -7r(r 2 2 +a 2 )~- 1 - 7r(r| + a 2 )^ 1 . (50) 

It is crucial to note that the action is independent of fa, and therefore we obtain the constrained 
instanton. The probability of the Kerr black hole creation is 

P k « cxp(7r(r 2 2 + a 2 )^- 1 + Tr(r 2 + a 2 )^ 1 ). (51) 
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It is the exponential of one quarter of the two horizon areas, or the total entropy of the universe. 



V. The Newman-de Sitter black hole 



Now let us turn to the charged black hole case. The vector potential can be written as 

Qr(dt - a sin 2 Odcf) +P cos 0(adt - (r 2 + a 2 )d(f>) 
A- — . (52) 

We shall not consider the dyonic case below. 

One can closely follow the neutral rotating case for calculating the action of the correspond- 
ing constrained gravitational instanton. The only difference is to add one more term due to the 
electromagnetic field to the action of volume. For the magnetic case, it is 

(53) 



(54) 



2S 2 V r'i + a 2 r\ + a* 
and for the electric case, it is 

h^iQ" 1 ( r 2 r 3 



2£ 2 \r 2 +a 2 r 2 + a 2 

In the magnetic case the vector potential determines the magnetic charge, which is the integral 
over the S 2 space sector. Putting all these contributions together one can find 

I = -Tr(r 2 + a 2 )-- 1 - 7r(r 2 + a 2 )-- 1 (55) 

and the probability of the creation of a magnetically charged black hole is 

P n « cxp(^(r 2 + a 2 )-" 1 + ir(r 2 3 + a 2 )-" 1 ). (56) 

In the electric case, one can only fix the integral 

(57) 



where the integral is around the S 1 direction. So, what one obtains in this way is ^f((j,a,hij). 
In order to get the wave function ^f(Q,a,hij) for a given electric charge, we have to repeat the 
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procedure like the Reissner-Nordstrom case. The Fourier transformation is equivalent to adding one 
more term to the action 



Then we obtain the same formula for the electrically charged rotating black hole creation as that 
for the magnetic one, 



It is easy to show that the probability is an exponentially decreasing function of the mass param- 
eter, charge magnitude and angular momentum, and the de Sitter spacetime is the most probable 
Lorcntzian evolution at the Planckian era. 

VI. Discussion 

The result of this paper has shown that the probability of the black hole creation is the exponential 
of the total entropy of the universe. The entropy is equal to one quarter of the sum of the black 
hole and cosmological horizon areas. 

The probability is an exponentially decreasing function in terms of the mass parameter, charge 
magnitude and angular momentum. Since this is only the confirmation of the conjecture, the result 
is no surprise. The only surprise is the fact that our result is independent of the choice of ji or f% 
for the formation of the constrained gravitational instantons. 

To get a meaningful result, one has to be careful to identify the meaning of the wave function; so 
for the rotating case and electrically charged black holes, one has to introduce Fourier transforma- 
tions into the calculation; otherwise the result becomes meaningless. It is interesting to note that 
Nature would give us a beautiful result if our request is reasonable. 

In quantum field theory, the temperature associated with a black hole is well defined. By using 
the reciprocal of the Hawking temperature as the period of the imaginary time, one can avoid the 
conical singularity at the horizon. However, if we remain only at thermodynamics level, and if 
one considers the reciprocal of the period for the constrained gravitational instanton as an effective 
temperature, then from the calculation, it seems the temperature can be taken quite arbitrarily. We 




(58) 



P n « exp(7r(rl + a 2 )E 1 + ir(r 2 + a 2 )Z r ). 



(59) 
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appear to overcome the obstacle that the temperature of the black hole and cosmological horizons, 
in general, are different. This makes our calculation feasible. Temperature is a very subtle concept 
even in special relativity, let alone in general relativity. A thorough discussion about temperature 
is beyond the scope of this paper. However, the concept of entropy is very clear in any case. 

Our calculation has also very clearly shown that the gravitational entropy is associated with 
topology of spacetime, as Hawking emphasized many times [18]. 

From the no-hair theorem, a stationary black hole in the de Sitter spacetime background can 
only have three parameters, mass, charge and angular momentum, so the problem of the quantum 
creation of a single black hole at the birth of the universe is completely resolved. 

Acknowledgment : 

I would like to thank Prof. G. Coyne of Specola Vaticana and Prof. R. Ruffini of Rome University 
for their hospitality. 

References: 

1. S.W. Hawking, Mon. Not. R. Astro. Soc. 152, 75 (1971). 

2. B. Carr and S.W. Hawking, Mon. Not. R. Astro. Soc. 168, 399 (1974). 

3. Z.C. Wu, Phys. Rev. D30, 286 (1984). 

4. L.Z. Fang and M. Li, Phys. Lett. B169, 28 (1986). 

5. Z.C. Wu, in Proceeding of the Fourth Marcel Grossman Meeting edited by R.Ruffini (North 
Holland, Amsterdam, 1986). 

6. R. Bousso and S.W. Hawking, Phys. Rev. D52, 5659 (1995). 

7. P. Ginsparg and M.J. Perry, Nucl. Phys. B222, 245 (1983). 

8. F. Mellor and I. Moss, Phys. Lett. B222, 361 (1989); Class. Quantum Grav. 6, 1379 (1989). 

9. I.J. Romans, Nucl. Phys. B383, 395 (1992). 

10. R.B. Mann and S.F. Ross, Phys. Rev. D52, 2254 (1995). 

11. S.W. Hawking and S.F. Ross, Phys. Rev. D52, 5865 (1995). 

12. J.B. Hartle and S.W. Hawking, Phys. Rev. D28, 2960 (1983). 

13. G. Hayward and J. Louko, Phys. Rev. D42, 4032 (1990). 

14. X.M. Hu and Z.C. Wu, Phys. Lett. B149, 87 (1984). 



15 



15. G.W. Gibbons and J.B. Hartle, Phys. Rev. D42, 2458 (1990). 

16. S.W. Hawking, Nucl. Phys. B239, 257 (1984). 

17. G.W. Gibbons and S.W. Hawking, Phys. Rev. D15, 2738 (1977). 

18. S.W. Hawking, in General Relativity: An Einstein Centenary Survey, eds. S.W. Hawking 
and W. Israel, (Cambridge University Press, 1979). 



16 



